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ABSTRACT 

We use N-body simulations to explore the origin and a plausible orbit for the Orphan Stream, 
one of the faintest substructures discovered so far in the outer halo of our Galaxy. We are able 
to reproduce its position, velocity and distance measurements by appealing to a single wrap of 
a double-component satellite galaxy. We find that the progenitor of the Orphan Stream could 
have been an object similar to today's Milky Way dwarfs, such as Carina, Draco, Leo II or 
Sculptor; and unlikely to be connected to Complex A or Ursa Major II. Our models suggest 
that such progenitors, if accreted on orbits with apocenters smaller than ^ 35 kpc, are likely 
to give rise to very low surface brightness streams, which may be hiding in the outer halo and 
remain largely undetected with current techniques. The systematic discovery of these ghostly 
substructures may well require wide field spectroscopic surveys of the Milky Way's outer 
stellar halo. 

Key words: galaxies: haloes - galaxies: formation - galaxies: evolution - galaxies: kinematics 
and dynamics. 



Tidal streams represent direct signatures of the merging history of 
galaxies. In Cold Dark Matter (CDM) models, structures grow in 
a bottom-up fashion, by the accretion of smaller sub-units. After 
entering in the gravitational domain of a larger system, tides effec- 
tively remove material from the satellites, creating "tails" of par- 
ticles, either dark matter, gas or stars that approximately trace the 
orbits of their progenitors. The existence of tidal streams in the 
halos of galaxies would therefore be a natural expectation in the 
hierarchical paradigm of structure formation. 

However, it is still unclear how fundamental or dominant 
mer gers have been in the build-up of our Galaxy. In a recent pa- 
per, Isell et al.l l20o3) found that the amount of substructure in the 
Milky Way's halo is consistent with the most extreme scenario, in 
which it was entirely formed from the accretion of satellites. 

Although a large amo unt of stellar streams are predicted in 



these hierarchical mo d els jjohnston et alj [T995 I: IjohnstonI 



Helmi & White 199 



gl: iBullock et alj I2OOII : Ikelmi et al 



199S 



2003 



Bullock & JohnstonI l2005 ir detection is generally difficult due 
to their low surface brightness and low contrast against Galac- 
tic field stars. This implies that our own Galaxy and also our 
closest neighbour, M31, may offer the best chances for identi- 
fying streams in stellar surveys. Perhaps the clearest and best 
known example is the Sagittarius Stream, discovered about a 



decade ago. With the adven t of large wide-field surveys like 

2MASS Is krutskie et al. 2006) and t he Sloan Sk y Digital Survey 

jYork et alJi200Q ; Adelman-McCarthv et al.ll2007i) it was possibl e 

to trace this stream out to 360° on the sky jMaiewski et alj|2004) , 

allowing the modelling of the orbit and intrinsic properties of its 

progeni tor, the Sagittarius dwa rf (Helmi & White 2001; Ibat a et al. 

1 200 ij : iMartinez-Delgado et al.ir2004l ; lHelmill2004l ; iFellhauer et al. 

l2006h . Also M31 shows a prominent "Giant Arc", a stellar 

stream of average surface brightness T,v ~ SOmag/arcsec'^, 

' — ' ' 1 

whose progenitor has not yet been confirmed (Ibata etal. 20011; 
Guhathakurta e t al. 2006 ; Kalirai et al. 2006; Fardal et al. 20oi; 
iFont et alj2006| - |Gilbert et al.l2007h . Therefore, these surveys have 
consolidated the idea that the stellar halos of M3 1 and the Galaxy 
may well be highly lumpy components. This appears to also be the 
case for galaxies beyond the Local Group, as recent studies have 
shown (e.g. iMartinez-Delgado et al..,2008i . and references therein). 



Stellar streams are useful tracers of their progenitor's orbit, 
due to their coherence in phase-space. This property allows us 
in some cases to link a given tidal feature with its parent satel- 
lite; for instance, by extrapolation of its great circle on the sky 
dLvnden-Bell & Lvnden-Belll 19951) . At the same time, tidal streams 
also hold essential clues about the object in which they originated. 
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Figure 1. Properties of the satellite model considered in this work, left: dark matter circular velocity (black squared dots, solid line) and the (projected) 
velocity dispersion of the luminous component (open red circles, dashed line) of the satellite as a function of radius. Lines indicate the profiles as set up in the 
initial conditions, while points show these quantities after the model is relaxed in isolation during 2 Gyr middle: dark and luminous matter density profiles for 
the satellite, arrows show the half-mass radii for each component. Lines and points are colour-codded as in the previous panel, right: projected positions of the 
dark matter (black) and luminous (red) pailicles in the relaxed satellite model. The green circles indicates the half-mass radius, highlighting the appreciable 
segregation of the stars with respect to the dark halo. 



as their luminosity and cross section are directl y related to th e 
mass and velocity dispersion of their progenitors ( |johnstonlll998h . 
Stream properties also depend on the time of disruption, since 
the material that is stripped off earlier in time tends to become 
broader, and hence, to ha ve lower surface brightness as time goes 
by (Helmi&Whitdl 19991) . In this context, we naturally expect an 
observational bias towards detecting remnants of recent accretion 
events involving massive progenitors, as the case of the Sagittarius 
stream in the Milky Way or the Giant Arc in M31. However, an 
unusually narrow and faint stream discovered in o ur Galaxy stands 
out from this general trend. The Orphan stream jGrillmaiil 1 200d : 
iBelokurov et"d]|2007l) with only ~ 0.70 kpc of full width half- 
mass (FWHM) projected on the sky, is about ~ 5 times narrower 
than the Sagittarius Stream, and has a surface brightness which is 
about a factor ~2 lower. Given its atypical properties, the Orphan 
Stream provides us with the opportunity of studying the fainter end 
of the family of objects that built up the stellar halo of our own 
Galaxy. 

As highlighted by its name, the Orphan Stream progenitor 
has not yet been identified, although its properties (mainly cross- 
section and luminosity) seem to favour a dwarf-like object rather 
than a globular cluster. Some attempts have been made to link 
this stream to other known objects of our Galaxy. Belokurov et 
al. (2007b), based on the path s defined by great circles on the sk y 
(a technique first developed bv lLvnden-Bell & Lvnden-Belllll99^ . 
highlighted a potential connection of the stream with several glob- 
ular clusters (Palomar 1, Arp2, Terzan 7 and Segue 1), as well as to 
the recently discov ered faint dwarf galaxy, Ursa Major II (UMa II, 
IZucker et al.ll2006h . In a later paper, Fellhauer et al. (2007) used 
numerical simulations to test the likely association between the Or- 
phan Stream and UMa II. These authors propose a model in which 
both objects have a common origin, also allowing for a physical as- 
sociation of the stream with a set of high velocity clouds named 
Complex A. However, this model requires the stream to be the 
result of the exact overlap on the sky of two independent wraps, 
which at face value appears somewhat contrived given its notice- 
able cohesion. 

In a novel approach, Ijin & Lvnden-Belll j2007h exploited the 
possible relation of Complex A with the stream to estimate its aver- 



age distance. However, the method under-predicts the distances by 
a factor of ~ 3, also disfavouring the scenario in which the Orphan 
Stream and Complex A share the same orbit. 

This paper presents a new attempt to find a plausible orbit for 
the Orphan stream progenitor, that is able to reproduce the current 
position and velocity measurements appealing to only one wrap; as 
suggested by the coherence observed in the images of the stream. 
Our model for the progenitor aims to be consistent with dwarf 
spheroidal galaxies, like those found orbiting the Milky Way today. 
The paper is organized as follows: we describe the simulations as 
well as the progenitor model in § |2l we present the results and dis- 
cussion in §[3]and §|4] Finally, we summarize our main conclusions 
in §111 



2 NUMERICAL MODELLING 

We use the N-body code GADGET-2 jSDringell2005l) nto study the 
evolution of a two-component satellite galaxy orbiting in the (fixed) 
gravitational potential of the Milky Way. We model the Galactic 
potential as follows: 

• a spherical NFW jNavarro et al.ll996l,ll997h dark matter halo: 



x{l + cxY 



with mass Mvir = 1 x 10^'^ M^, concen tration c = 12 and virial 
radius r^ir = 258 kpc jKlvpin et al.ll2002h . 
• a Hemquist bulge: 



p(r) = 



27rr (r + at,)'^ 



with mass Mug = 3.4 x 10 Mq and scale length at — 0.7 kpc, 
• a Miyamoto— Nagai disk: 



^ http://www.mpa-garching.mpg.de/~volker/gadget/index.html 
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p{R,z) = 



27r [R-2 + (a + + 62)2)5/2(^2 + ^2)3/2 



with parameters: Afdsfe = 1 x 10^^ M©, a = 6.5 kpc and b — 
0.26 kpc (Johnston et al. 1999). 

The composite circular velocity at the solar distance is Vc = 
220 km/s in agreement with observations. The circular velocity at 
the virial radius is Vc(r„i,.) = Vvir ~ 136 km/s. Note that for sim- 
plicity we have assumed a spherical dark matter halo, and that the 
mass distribution is time-independent. This assumption is justified 
because the gravitational potential inside the orbit of the Orphan 
stream (< 40 kpc, see next Section) is unlikely to have changed sig- 
nificantly in the recent past becau se of the presenc e of a r elatively 
old thin disk. On the other hand IPenarrubia et al. I ( l2006h suggest 
that the properties of tidal streams mainly reflect the present-day 
potential of the primary halo and are not fundamentally affected by 
its growth in time. 

Our model of the satellite has two spherical components: an 
extended dark matter halo and a more concentrated "luminous" 
component (see middle and right panels of Figure [T](. Dark mat- 
ter particles are distributed following a Hemquist profile of mass 
MJ^fe = 2.5 X 1O*M0 and characteristic scale asat = 0.9 kpc. 
The circular velocity of the dark halo peaks at Vmax ~ 0.75 kpc 
reaching Vmax ~ 20 km/s (see left panel of Figure [T), and hence 
would be consistent with the properties of the Milky Way dwarf 
spheroidals (Pefiarrubia et al. 2007,2008). The second component 
mimics the stellar content of a dwarf galaxy and is represented by a 
Plummer profile of mass M^^^ — 7.5 x 10^ M0 and characteristic 
scale bsat =0.1 kpc. The half light radius of our model satellite is 
rh ~ 0.13 kpc and its central velocity dispersion is ao ~ 7 km/s. 

We used the web-tool BaSTl[3to generate isochrones and con- 
vert the stellar mass {Mstr ~ 7.5 x 10'' Af©) to the luminosity 
of our modelled satellite. Given the typically old stars present in 
dwarfs galaxies, and assuming a single population of 10 Gyr and 
[Fe/H]~ —2 dex, the conversion factor we obtain is 7 ~ 2.9. 
This gives a total luminosity for the progenitor in our simula- 
tions of L ~ 2 X 10^1/0, consistent with the luminosity of dwarf 
spheroidals in the Local Group. This object is dark matter domi- 
nated, with a mass-to-light ratio 7 ~ 39 measured within the half 
light radius. With these parameters the model is comparable to sev- 
eral classical dwarfs, such as Carina, Leo II or Sculptor. 

The number of particles used in our simulations are 5x10''' 
and 2x10^ for the dark and stellar components, respectively. We 
use an unequal softening scheme, with a Plummer-equivalent soft- 
ening length: eoM = 0.04 kpc and estr ~ 0.004 kpc for the dark 
matter and luminous component, respectively. We set up th e initia l 
conditions following the procedure outlined bv lHemauisJjl993h . 
aimed to generate multi-component systems in dynamical equilib- 
rium. The satellite is first evolved in isolation during 2 Gyr (~ 200 
crossing times for a particle at the half-light radius, while for one 
at the edge, i.e. at 5 kpc from the center, this corresponds to ~ 5 
crossing times), where it is allowed to relax. Then, the satellite was 
put on its current orbit in the fixed Milky Way potential. 

As stated above, our model for the Orphan stream pro- 
genitor could well represent one of the "classical" Milky Way 
dwarf galaxies, where "classical" is used to distin guish the first 
eleven discovered s atellites around our Galaxy jMateol 1 19981 : 
Ivan den Bergdll999h . from the more recently identified SDSS 

2 www.te.astro.it/BASTI/index.php 
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Figure 2. Distribution of dark matter (black) and star particles (red) at 5.3 
Gyr after infall. The section matching the Orphan Stream measurements 
is enclosed by the green rectangle. It sits on the trailing arm, traced here 
(together with the leading arm) by the dashed blue line. The cyan asterisk 
corresponds to the position of the most bound star particle, with the sense of 
motion in the orbit shown by the arrow. Dark matter particles are removed 
from the progenitor at earlier times, causing the stream to be dynamically 
old and with poor coherence in space, although characteristic "shells" or 
"caustics" can still be distinguished in the distribution. On the other hand, 
star particles are more clumped, approximately tracing the last loop {trailing 
arm) and next loop (leading arm) of the orbit. 



dwarfs iZucker et al .120041 : 1 Willman et alj2005l:lMartin et al. | 2o6i 



Belokurov e t al. 2006* , Zucker et al 2006a b iBelokurov et alj2007 : 
Irwin et alj ,2007 : .Maiewski et alj |2007| : Isimon & Gehal l2007h . 
These new dwarfs are typically ~ 100 times less luminous, but with 
comparable half light radii (~ 80 — 500 pc) and velocity disper- 
sion (~ 4 — 8 km/s). Therefore their surface brightness is apprecia- 
bly lower (~ 28 — 30 mag/arcsec'^) and their mass-to-light ratio is 
~ 10 — 50 times higher. One of the strongest constraints on the Or- 
phan Stream progenitor comes from its total luminosity. Belokurov 
et al. (2007b) estimate that the total r-band absolute magnitude for 
the (visible) portion of the stream is Mr ~ —6.7. Since leading 
and trailing arms should have comparable mass, a lower limit to 
the total luminosity is roughly L ~ 10^ Lq, ruling out practically 
all new SDSS dwarfs as suitable candidates 0. 



3 RESULTS 

It is not possible to fully constrain the orbit of the Orphan stream 
with only the available radial velocity and distance information. 
Proper motion measurements would be needed to that end. How- 
ever, given that the stream extends ~ 55° on the sky, it is possible 



^ Only the Canes Venatici I dwarf, whose luminosity L ~ 1.2 x 10^ L0 
is above this lower limit. However, this galaxy is too extended (r^ = 0.56 
kpc) to give rise to a stream as narrow as observed 
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to find a suitable orbit by requiring that it should pass through both 
ends of the observed stream, with the measured (although prelim- 
inary) radial velocities and distances. This can be done by requir- 
ing that the total angular momentum of both ends of the stream be 
the same (for example the orientation of the orbital plane is con- 
strained by the cross product of the position vectors of the stream 
end-points). With this condition, we randomly generate possible 
orbits that match (within the errors) all observables. 

The satellite is placed at one apocenter of such an orbit, from 
where we follow its subsequent evolution within the Milky Way's 
gravitational field. The orbit is confined to the inner 40 kpc of the 
Galaxy halo (apocenter: Vapo ~ 38 kpc, pericenter: rper ~ 7 kpc); 
where the gradient of the potential is large, producing strong tides 
that fully disrupt the progenitor in less than ~ 3.5 Gyr. The snap- 
shot view with the final distribution of satellite particles after 5.3 
Gyr is shown in Figure |2] The box is 100 kpc on a side and the 
projection is along the x-coordinate. The green rectangle shows the 
portion of the stellar trail that would represent the Orphan stream in 
our model. The distribution of dark matter (black dots) is broad and 
shows no spatial coherence. Due to its more extended distribution, 
dark matter particles are stripped off first, and hence this debris is 
dynamically older, and consequently more diffuse. However, we 
can still identify by eye a few she lls or "caustics" that outline the 
turn-around points alo ng the orbit jHemquist & Ouiririlll987l , ri988l : 
iHelmi&WhiSi 1999b 

The luminous component of our satellite is a factor ~ 17 more 
concentrated (the ratio of their half mass radii is: r'f,^^ /r^^ — 
16.7), which makes its core more resilient to tidal disruption than 
the extended dark matter. Nevertheless, after the second pericenter 
passage, the first stellar trails are produced (at which point the satel- 
lite has lost more that 95% of its initial dark matter content). The 
stellar streams hence grow at the expense of the satellite's pruning, 
a process that, if strong and enduring enough, drives to the irre- 
versible disruption of the object. 

Note that we have followed the evolution of the model satellite 
for 5.3 Gyr. This timescale is partly driven by the initial properties 
of the satellite and by our desire to reproduce the current proper- 
ties of the stream. As discussed in detail in § 13.31 this choice is 
not unique as strong dependencies between the initial conditions 
of the model satellite and the integration time exist. Nonetheless, 
the above choice of an integration timescale shorter than a Hub- 
ble time for a relatively tight orbit may be also justified in light of 
recent cosmological simulati ons. These have sh own that satellites 
are often accreted in groups d Li & Helmi|[2008l ). and that through 
multiple-body interactions acting during the tidal diss ociation of 
such groups their orbits can be changed significantly jSales et al] 
l2007l : lLudlow et alj2008h . Such events are relatively common, and 
might explain how our model satellite was put on its current orbit 
only 5.3 Gyr ago. 

3.1 The orbit and the stream global properties 

Figure [3] shows the position on the sky in galactic coordinates, 
the heliocentric distances and velocities of the star particles in our 
model, in comparison with the measurements of the Orphan Stream 
reported by Belokurov et al. (2007b). 

In this snapshot, the former center of mass of the satellite is 
located at 1 = 306.1°, b=13.9°, at a heliocentric distance of 13 kpc 
and moving with a velocity of vq = 97.6 km/s. 

The most recent wrap of the orbit is enclosed between galac- 
tic latitudes b=[-30°,+60°], and is pretty well traced by the majority 
of the star particles (see upper panel in Figure[3]l. However, we can 
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Figure 3. All sky view of the final distribution of stellar particles in our 
model. We show in the upper panel the projected positions in galactic co- 
ordinates {l,b), middle and bottom panels correspond to the heliocentric 
distances {Dq) and velocities (Vq), respectively. The thin red line shows 
(a portion of) the trajectory of a point-mass on the same orbit as the Oiphan 
Stream progenitor, integrated in the analytic potential of the Milky Way. 
Green crosses indicate the measurements for the Orphan Stream, fields 1-5 
taken from Belokurov et al. (2007). For comparison with previous works, 
we also include current estimations for the positions, distanc es and veloc- 
ities for Complex A (blue e mpty circles) and Ursa Major II iMartin et alj 
l200l [s imon & Gehj|2007l magenta filled triangle,). Red asterisks show 
the properties of 5 red giant stars from the Spaghetti survey, probably re- 
lated to the stream (see §|4]for further details). 

also see evidence of previous loops, visible as a bulk of particles 
with 1< 90° and also at 1> 280° . The section of the stream that 
matches the Orphan Stream observations forms an arc of ~ 55° 
in the trailing arm (the sense of motion in the orbit is pro-grade, 
i.e., towards larger galactic longitudes), where the maximum sur- 
face brightness is reached. The leading arm is responsible for the 
overdensity of particles with 1< 90° and also 1~ 300°. The differ- 
ent branches are better distinguished in the middle panel of Figure 
[3] where the leading arm has a maximum distance of ~ 40 kpc and 
the particles in the tip of the trailing arm are located at 1> 320° 
with distances up to Dq ~ 50 kpc. 

Figure [3] shows that positions, distances and velocities de- 
termined for fields 1-5 (Belokurov et al. (2007b)) of the Orphan 
Stream are well reproduced by our model. Our distances might be 
slightly larger than measured by Belokurov et al. (2007b) in fields 
2 and 3, but still in good agreement, especially given the size of 
the error bars. The velocities are also nicely consistent with the 
observed values, considering that the latter were determined with 
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Figure 4. Transverse surface brightness profile of the simulated stream, in 
the region that matches the location of the Orphan Stream. For each longi- 
tude range we re-centred the distribution according to its baricenter latitude, 
f)0 ■ The main box shows the total distribution of light, while the small in- 
dicates the likely "detection", i.e., the excess of flux over a background 
surface brightness = 34 mag/arcscc^. As discussed in the text, ac- 
counting for the foreground and background contribution helps to reproduce 
the narrow cross-section of the Orphan Stream shown in SDSS data. 



barely more than a couple of dozen stars per bin. Although some 
attempts to improve the observational estimations of the velocities 
in the Orphan Stream fields have already been made, unfortunately 
they have not yet succeeded (Belokurov, private communication). 
Nevertheless, it is clear that a more definitive test to the orbit (and 
the progenitor) presented in this work will come from better mea- 
surements of velocities and smaller distance brackets for each of 
the fields. 

A possible association between the Orphan S t ream and Com- 

S ilex A was suggested by Belokuro v et al.l ( |2007|) : iFellhauer et al.l 
20071) : IJin & Lvnden-Belll l l2007il . Complex A is a set of high 
velocity clouds whose position in the sky lies along the best- 
fitting great circle trac ed by the Orphan Stream (see Figure 7 in 
iBelokurov etai] ( l2007l) ). However, the model presented here shows 
no connection (in the current as well as in previous loops) between 
these two objects (see Fig.[3]l. Our results also do not support an 
association between the stream and UMa II, as suggested by Fell- 
hauer et al. (2007). 



3.2 Stream intrinsic properties 

In the following paragraphs we address how well the intrinsic prop- 
erties of the Orphan Stream (magnitude, extension, width, etc) are 
reproduced in our model: 



• Stream width 

One of the peculiarities of the Orphan Stream is its small cross sec- 
tion (~ 2°), that distinguishes it among its siblings, e.g., the Sagit- 
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Figure 5. Main box: All-sky surface brightness map of the simulated 
satellite and its debris, assuming a mass-to-light ratio for star particles of 
7 = 2.9 (see text for details). Small box: zoom-in into the Orphan Stream 
region. Notice the change of scale in the density map. In this box we have 
considered = 34 mag/ arcsec^ as the lower limit in the surface bright- 
ness, in order to account for the foreground and background stars on the 
Milky Way in this region of the sky. The average FWHM of the "observ- 
able" portion of the simulated stream is then ~ 2° in nice agreement with 
the observations. The trend of increasing number of stars towards larger 
latitudes is also reproduced in our model. 



tarius Stream or Monoceros Ring. The total width of the model 
stream in the region 1=[200° , 250°] is 15°, significantly wider than 
what has been reported. However, much of this stream is at such 
low surface brightness levels that it would remain undetected given 
the expected background, as shown for example in Figure 5. To ac- 
count for the contribution from background and f oreground stars 
we use the Besancon ModeJ3 4 l lRobin et alj|2003b . From this we 
estimate the number of stars in our Galaxy that are expected within 
a 15 square degree region of the sky centred at the Orphan Stream's 
position. We apply the apparent magnitude limit and color cuts 
adopted by Belokurov et al. (2007b) (stars with 21 < r < 22 
and — r) < 0.4): obtaining an r-band surface brightness for the 
field: Ef ~ 33.8 mag/arcsec^. Therefore, the stream stars would 
be confused with the background at Er ~ 34mag/arcsec^. When 
we apply this threshold to our simulation, we find that the "ob- 
servable" portion of the stream nicely reproduces the ~ 2° cross- 
section. This is illustrated in detail in Figure 4, where we show the 
average cross-section profiles of the simulated stream and the effect 
of adding the foreground and background contributions. The final 
structure of the stream is shown in the bottom panel of Figure 5. 



http://bison.obs-besancon.fr/modele/ 
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• Stream length 

The current detection of the Orphan stream covers an arc of ~ 55° 
on the sky. This strongly constrains the time required for the debris 
to spread in the simulation to a matching extent. For our model we 
need about ~ 5 Gyr. Shorter timespans of the modelled satellite 
in the fixed potential {t < 4.5 Gyr) generate streams that are not 
long enough, while a longer time integration (t > 6.3 Gyr) creates 
tidal streams that are too wide and diffuse compared to the Orphan 
Stream. 

• Stream luminosity 

The total amount of mass in the simulated stream depends on the 
limit in the surface brightness we adopt. Applying the threshold 
given in the observations by the field brightness, S^, we obtain 
m ~ 1.1 X W^Mq. This translates to a luminosity: L ~ 2.3 x 
10*1/0, or an absolute magnitude AU ~ —6.4. These numbers are 
in very good agreement with Belokurov et al.'s measurements, who 
estimate Mr ~ —6.7. Nevertheless, we notice that the real mass in 
the stream is twice this number, once the limiting surface brightness 
considered drops to Er = 38mag/arcsec^ instead of the = 
34mag/arcsec^ imposed by background and foreground stars. 

• Stream gradient 

Our model also reproduces the non-uniformity in the surface 
brightness of the Orphan Stream. As pointed out by Belokurov 
et al., the density of stars noticeably increases towards lower lati- 
tudes. Bottom panel of Figure|5]shows that there are approximately 
twice as many stars in 1 > 230° (fields 1-2) than for the upper 
tip 1 < 220° (fields 4-5). This effect is generated by: i) fields 1-2 
are closer to the former progenitor than the upper tip of the stream, 
and ii) fields 4-5 are 50% farther away from the Sun than the lower 
section of the stream, compromising even more its detection. 

• Final fate of the progenitor 

The large degree of disruption in our model satellite needed to gen- 
erate an arc on the sky as large as ~ 55° prevents the progenitor to 
survive as a self-bound entity for more than f ~ 3 Gyr. A robust 
conclusion that we can draw after having explored several mod- 
els for the satellite, starting from conditions that resemble today's 
Milky Way dwarfs, is that any remnant at the present time should 
be very close to (if not already) fully disrupted. 

3.3 Model degeneracies 

The orbit presented here is likely to be one of many possible ones. 
This is because of the freedom introduced by the lack of proper 
motions measurements and the large uncertainties in the radial ve- 
locities and distances, as well as by the unknown exact form of the 
gravitational potential of the Galaxy. 

Furthermore, as a final comment in this section, we point out 
that the uncertainties on the measurements and the lack of an entire 
sky survey with the resolution of the SDSS, makes the problem of 
determining the properties of the Orphan Stream progenitor quite 
degenerate (see Fellhauer et al. (2007) for a similar conclusion). 
The approach in this paper has been to choose a object that resem- 
bles one of today's Milky Way dwarfs. Nevertheless, the solution 
is not unique. 

There are two factors that largely determine the success of a 
model on the orbit presented here: the central density of the satellite 
(or the average density within the half-light radius) and the scale of 
the stellar distribution (see for instance Figure 2 of Penarrubia et al. 
2008). Note that these two factors also determine the stellar velocity 
dispersion, a quantity that has direct impact on the stream's width. 

The satellite's central density determines the likelihood of sur- 
vival and the degree of disruption of the object on the current orbit. 



Hence, this will determine how long the streams are and for how 
long the object has been orbiting the Galaxy. For example, the ob- 
servations constrain the length of the stream to an arc of 55 deg. 
Assuming that the Orphan Stream is actually the brightest section 
of the total (full-sky) stream, and given that no progenitor has been 
found close to the region where the stream reaches its peak surface 
brightness, this implies that the object must be close to (if not) fully 
disrupted. Therefore the initial average density of the object must 
be comparable to the density of the host at the pericenter distance 
(~ 3 X 10^ A/© kpc^'^). An object with a significantly higher cen- 
tral density will not suffer enough tidal stripping nor shocks to be 
significantly perturbed and disrupted. An object with a much lower 
density will be fully disrupted giving rise to streams that are too 
faint to be observable. This is why a model like ours needs 5.3 Gyr 
to give rise to an object that is fully disrupted and to streams of the 
required length. A model with a slightly higher dark-matter density, 
but the same stellar distribution, would survive longer, and hence 
the integration time required to reach the present state would also 
need to be larger. A model with a smaller amount of dark matter 
near the center (e.g reduced by 50%) is also possible, since it does 
not lead to significant changes in the observables, as long as the 
satellite evolves for a shorter time in the Milky Way potential (4.7 
Gyr instead of 5.3 Gyr). 

The second factor that is crucial in our model is the extent of 
the luminous component, since this determines almost completely 
the narrowness of the observed streams. In principle, an interplay 
between a shorter time integration and a more extended progenitor 
could also provide a suitable model for the stream. However, we 
were unable to obtain a good match between simulations and obser- 
vations by using a satellite model as extended as the one proposed 
by Fellhauer et al. (2007). These authors used a two-component 
model with similar mass for the luminous particles, but its half light 
radius is ~ 480 pc, approximately 3.5 times larger than our fidu- 
cial progenitor. Such an extended object in our preferred orbit pro- 
duces a stream that is too wide and with too low surface brightness 
in comparison to the Orphan Stream. On the other hand, the mass 
in the luminous component could also be increased maintaining 
the mass-size relationship observed for the local dwarf spheroidals. 
However, these more extended progenitors give rise to tails that are 
too wide and (now) too bright to fit the Orphan Stream observa- 
tions. 

As stated above, these degeneracies are important since the 
integration timescale can then be made as long as desired by in- 
creasing the average density of the progenitor. For example, one 
may desire a model which has been on the present orbit around the 
Galaxy for the past 10 Gyr. By performing numerical experiments 
we have found that a satellite with a mean density within the half 
light radius which is 10 times larger than our fiducial model 
(i.e., ~ 2.9 X IO^Mq kpc"^) will achieve after 10 Gyr the level of 
disruption needed to reproduce the observations. Such a model can 
be obtained by increasing the dark matter halo mass by a factor 10 
while keeping its scale-length asat fixed. However, this will lead 
to a significantly larger stellar velocity dispersion in comparison to 
the previous model (with ~ 7 km/s), and hence will fail in repro- 
ducing the width of the Stream. Therefore also the extent of the lu- 
minous component must be modified to become more concentrated 
than our fiducial model. The magnitude of this modification can 
be derived from the following argument based on the initial stellar 
velocity dispersion of the system astr- Since most of the mass is 
in dark matter, we can neglect the contribution of the stars to the 
velocity dispersion, and therefore assuming, as before, a Hemquist 
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halo and a Plummer profile for the stars we find: 

where Uaat and bsat are the charact eristic radii of th e Hemquist and 
the Plummer profiles, respectively i Hemquis^l 1993h . The dark halo 
mass, Mdrk is fixed by our requirement of a given central density. 
If we now vary bsat we can get a range of velocity dispersions 
that goes from 4 km if bsat ~ 20 pc up to ~ 24 km s^^ 
for bsat =0.1 kpc as before. If we use our fiducial model as a 
guide, the new model integrated for 10 Gyr will be successful if 
o'str ^ 7 km s^^ after ^ 5 Gyr of evolution. Since disruption 
causes the velocity dispersion to drop, models with initial astr in 
the range 9~I2 km s^^ will generate the kind of progenitor we are 
looking for. According to the equation above, that corresponds to a 
half light radius th ~ 50-65 pc. 

This new progenitor should also fulfil the constraints in lumi- 
nosity given by the brightness of the visible portion of the Orphan 
Stream (see Section 2), which implies L 10^1/0. Therefore, we 
find that by fixing the central density, and changing the parameters 
of the luminous component accordingly, we can generate a family 
of progenitors that are roughly consistent with today dwarfs satel- 
lites. On this basis, we conclude that models evolving for 5 to 10 
Gyr and with the properties specified above could represent feasible 
candidates to the progenitor of the Orphan Stream. 

For completeness, we have also explored single-component 
(no dark matter) models (reminiscent of a globular cluster or a 
dwarf galaxy that has already lost its subhalo) as tentative pro- 
genitors for the stream. We find that such objects should be more 
concentrated (typically half light radius ~ 40 pc) and need to be 
integrated for a longer time (~ 6.5 Gyr) than our fiducial model, in 
order to reproduce the length and width of the Orphan Stream. As 
discussed above, a higher density progenitor could be integrated 
for a longer timescale leading to a similar configuration. It is in- 
teresting to note that an object like a globular cluster has too high 
a central density (~ 10^^i\fo kpc""^) to be fully disrupted within 
a Hubble time on this orbit. The amount of mass in the currently 
detected stream is fixed from observations, hence, the total lumi- 
nosity of these models should also be L > 2 x 10^ Lq, as in the 
double-component case. Progenitors with such properties could be 
considered unlikely given the properties of the satellites in the Lo- 
cal Group: they are too bright and extended to be a globular cluster 
and too compact to represent a Milky Way dwarf galaxy (see for 
instance Fig. 10 of Peiiarrubia et al. 2008). Therefore, although we 
cannot rule out a single-component progenitor with such character- 
istics (since today's satellites might not be representatives of ear- 
lier accreted ones); the lack of observations of similar objects turns 
the dark-matter dominated model presented before a more suitable 
candidate for our analysis. 



4 PRELIMINARY PREDICTIONS AND FUTURE 
DETECTION PROSPECTS 

The upper panel of Figure |5] shows a dense enough section of the 
stream, likely identifiable with the current sensitivity of the SDSS; 
although unfortunately, sitting outside its coverage area. This new 
stretch contains about ~2 times less stars and traces an arc of 
~ 35° following the current Orphan Stream detection. It extends 
from 1 = 250° to 1 ~ 290° in the North galactic hemisphere. 
The leading arm could in principle be also identified in a survey 
covering the southern galactic hemisphere. However, its position 



towards the galactic center (1 = [—60°, 70°]) and its lower lati- 
tude (b — [—30°, 10°]) might compromise the chances of positive 
imaging detection. In addition to this, its mean surface brightness 
is on average lower than that of the trailing arm, where the Orphan 
Stream sits. This difference can be attributed to the fact that the por- 
tion of the trailing arm identified as the Orphan Stream lies close 
to apocenter of the orbit, where the spacial coherence of particles 
is locally enhanced (e.g., formation of shells or caustics seen in 
Figure |2} reaching the maximum surface brightness across the full 
sky. 

An interesting alternative to these photometric approaches, is 
offered by studies involving also stellar kinematics. Stars remain 
coherent in velocity space as a consequence of the conservation of 
phase-space density (Helmi & White^l999.) . enhancing the correla- 
tions of their dynamical parameters. 

One relevant study of stellar kinematics of (outer) halo stars 
is provided by the Spaghetti Survey. This pencil-beam survey of 
high-latitude fields, first described in Morrison et al. (2000), uses 
solely red giants as distant halo tracers. The final dataset consists 
of 102 spectroscopically confirmed giants (Starkenburg et al. 2008, 
in preparation). We find five giants within this dataset that could 
be related to the stream in our model (shown by the red asterisks 
in Figure[3](. Three of them are located right on top of the Orphan 
Stream fields as detected in Sloan. For two of these, the projected 
positions, distances and also heliocentric velocities in good agree- 
ment with the simulations, indicating a large probability of a real 
physical link to the stream. The case for the third giant is less clear, 
since its position agrees with the model , but its velocity is too l arge 
(given the radial velocities measured bv lBelokurov et all ( |2007|) are 
still preliminary, it is possible that this giant is physicEiIIy related 
to the Orphan Stream, but not in our model). The remaining two 
giant stars could belong to different wraps, although they sit on 
relatively low-density portion of the stream. Table[T]quotes the po- 
sitions, velocities and metallicities of the selected giants measured 
by the Spaghetti survey, together with our assessment of the likeli- 
hood of membership to the Orphan Stream progenitor according to 
our simulations. 

While red giant surveys like the Spaghetti Survey may be suit- 
able to find stream members and constrain the stream's proper- 
ties such as radial velocity, the number of detectable candidates 
is strongly limited by the surface brightness of the stream. The ap- 
proximate surface brightness of ~ 32.4 mag arcsec""^ (Belokurov 
et al. 2007b) tr ansforms to barely 1.3 giants per square degree 
( lMorrisorJl993h . With more fields on the stream and a more careful 
analysis of the existing SDSS spectra for giant candidates in fields 
on the stream, it should be possible to constrain the radial velocities 
better. However, even fainter substructures might need other tracers 
such as main sequence stars to determine their properties. 

Note also that some sections of the stream in our simu- 
lations fall relatively close to the sun (heliocentric distances of 
the ord er of 5 kpc) and co uld, in principle, be accessible to 
RAVE jSteirmietz et alj|2006l) . Measurements of the proper mo- 
tions of the stars associated to the stream could also be used 
to falsify our model. The orbit presented here as {fia,lJ-s) ~ 
(-0.10, -2.42)mas/yr at q = 162° and S = 0°, i.e. on Field I. 
These are within reach of future astrometric missions like GAIA. 

5 SUMMARY AND CONCLUSIONS 

We have presented N-body simulations that reproduce the prop- 
erties of the recently discovered Orphan Stream. We have stud- 
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Table 1. Galactic coordinates, distances, velocities and metallicities of the Spaghetti Survey giants shown in Figure|3] Errors in the quantities are also quoted. 
The last column denotes how consistent the membership to the Orphan Stream progenitor is according to our model. 



Location I b Dq ^0 '^^o Membership 

(deg) (deg) (kpc) (kpc) (km/s) (km/s) (dex) (dex) 



Field 3 


223.45 


43.44 


29.71 


3.71 


71.8 


5.5 


-1.38 


0.27 


likely 


Field 2 


233.83 


53.71 


23.46 


2.03 


9.6 


7.5 


-1.24 


0.26 


likely 


Field 1 


252.28 


52.98 


16.72 


1.39 


137.1 


10.1 


-1.24 


0.26 


unlikely 


Older wrap 


99.38 


-47.59 


33.88 


7.37 


-129.9 


14. 


-1.60 


0.59 


possible 


Older wrap 


333.34 


46.50 


37.16 


5.10 


168.7 


10.1 


-1.72 


0.27 


possible 



led an orbit able to match, the positions, velocities and distances 
of the stream. The orbit is in good agreement with the continu- 
ous nature of the stream, which may well be a single wrap of the 
trailing arm of a fully disrupted satellite. The main features of the 
stream, such as length, cross-section, luminosity and surface bright- 
ness are also recreated successfully in the model. However, the 
orbit does not provide an obvious association of the stream with 
Complex A or U rsa Major II, as had been suggested previously by 
Ijin & LyndeiTBelli ( |2007|) or Fellhauer et al. (2007). 

The satellite in our model consists of two components: an ex- 
tended dark matter halo within which the stars are deeply embed- 
ded. We find that the Orphan Stream progenitor could have been 
an object that resembles one of today's Milky Way dwarfs, such 
as Carina, Leo II, Draco or Sculptor, and hence, likely to be dark 
matter dominated. The model predicts the total disruption of this 
object in order to generate the ~ 55° stellar stream needed to 
match the observations. This indirectly implies that objects with 
the properties of the above-mentioned dwarfs can only survive as 
self-gravitating systems today, because they populate the outer re- 
gions of the Galaxy potential, where the gravitational forces are 
smaller, and the dynamical times of the orbits are longer. 

The Orphan Stream distinguishes itself from the previously 
known tidal features (e.g., Sagittarius Stream or Monoceros Ring 
in our Galaxy, or the Giant Arc in Andromeda) by its small cross- 
section. With an average width of ~ 2° it is about ~ 5 times nar- 
rower than the Sagittarius Stream. However, as argued above, this 
does not imply that its progenitor is an exceptional object. Such 
streams may therefore be far more common than suspected thus 
far. The visible strip of the stream may also be just the "tip of the 
iceberg", where the amount of mass under our detection limit could 
be at least twice the current estimation. 

We conclude that additional discoveries of tidal features such 
as the Orphan Stream are fundamental to our understanding of the 
outer halo of galaxies, as well as to tests of current cosmological 
models. The detectability of such coherent but very faint struc- 
tures is compromised not only by the initial luminosity of their pro- 
genitors, but also by the dynamical ages of these structures, since 
they become fainter as they age and diffuse away. The design of 
techniques to recover these ghostly features needs to be carefully 
thought out, in order to enable the detection of streams from ob- 
jects originally populating the faint end of the luminosity function. 
Future panoramic surveys of the Milky Way stars, as well as the 
determination of their kinematics, ages and metallicities will play 
a fundamental role in the Galactic astronomy of the next decades, 
helping to unveil the secrets recorded in the stellar halo of our own 
Galaxy. 
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